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COMPARISONS OF TWO-DIMENSIONAL 
SHOCK-EXPANSION THEORY WITH EXPERIMENTAL 
AERODYNAMIC DATA FOR DELTA-PLANFORM WINGS 
AT HIGH SUPERSONIC SPEEDS 

By Lloyd S. Jernell 
Langley Research Center 

SUMMARY 

An investigation has been conducted to explore the potential for optimizing airfoil 
shape at high supersonic speeds by utilizing the two-dimensional shock-expansion method. 
Theoretical and experimental force and moment coefficients are compared for four delta - 
planform semispan wings having a leading-edge sweep angle of 65° and incorporating 
modified diamond airfoils with a thickness-chord ratio of 0. 06. The wings differ only in 
airfoil maximum-thickness position and camber. The experimental data are obtained at 
Mach numbers of 3. 95 and 4. 63 and at a Reynolds number of 9. 84 x 10^ per meter. 

Although maximum lift-drag ratio is overestimated by approximately 6 to 8 percent, 
theory provides an accurate estimate of the relative effects of airfoil maximum -thickness 
position and camber. A relatively simple method is developed for predicting, in terms 
of lift -drag ratio, the optimum modified diamond airfoil at high supersonic and hypersonic 
speeds. 


INTRODUCTION 

Two-dimensional shock -expans ion theory has been widely used to determine the 
aerodynamic characteristics of sharp-leading-edge wings. References 1 to 3 present the 
results of analytical investigations which have employed the shock- expans ion method to 
determine airfoil characteristics at high supersonic and hypersonic speeds. The data of 
reference 3 indicate that this method provides an accurate estimate of the surface pres- 
sure up to arbitrarily high Mach numbers provided the flow acts as an ideal gas and the 
flow deflection angles are somewhat less than those initiating shock detachment. With 
regard to a wing of finite span at high Mach numbers, exact numerical solutions to the 
problem of flow over a flat delta wing with supersonic leading edges are developed for the 
upper and lower surfaces in references 4 and 5, respectively. In both investigations it is 
concluded that the normal force is practically independent of sweep angle and may be 
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approximated by two-dimensional theory. In reference 6 experimental force and pressure 
data for flat -bottom wings having wedge airfoil sections and sweep angles of 70° and 76° 
are compared with various theories at a free-stream Mach number of 8. 2. Again, the 
data indicate that the force coefficients may be closely approximated by the two- 
dimensional shock-expansion theory. 

The effectiveness of the shock-expansion method and other theories in predicting 
the pressure distribution over a wing-body configuration at high supersonic speeds is 
reported in reference 7. The wing is of delta planform with a sweep angle of 65° and 
employs a symmetrical diamond airfoil. The data show that at a Mach number of 4.63 
and moderate angles of attack (to approximately 11°), the experimental pressure distribu- 
tions over the wing upper surfaces are practically constant for each flat surface and agree 
well with the predictions based on two-dimensional Prandtl-Meyer expansion. Although 
considerable pressure gradient exists over the lower surfaces, the average pressure 
coefficient for each flat surface differs from the two-dimensional shock-expansion pre- 
dictions by less than 10 percent. In addition, the data show that as Mach number is 
increased, there is a decrease in the pressure gradient and an improvement in the agree- 
ment between theory and experiment. 

The foregoing comparisons suggest the possibility of employing the shock -expans ion 
method to optimize airfoil shape for a finite wing at high supersonic and hypersonic 
speeds. The primary purpose of this investigation is to explore the potential for such 
optimization in terms of lift-drag ratio in the high supersonic speed range. The study is 
limited to a series of modified diamond airfoils in conjunction with delta-planform wings 
of 65° leading-edge sweep. A relatively simple method is subsequently developed for 
predicting the optimum modified diamond airfoil, according to the two-dimensional shock- 
expansion theory, at high supersonic and hypersonic speeds. Theoretical predictions are 
compared with experimental force and moment coefficients obtained on several semispan 
wing models which represent a range of airfoil geometry variables spanning the optimum 
profile. The optimum wing is designed for a free-stream Mach number of 4.63. Experi- 
mental data are obtained at Mach numbers of 3, 95 and 4. 63 and at a Reynolds number of 
9. 84 x 10® per meter. 


SYMBOLS 


c chord 

c wing mean geometric chord 

Op drag coefficient, Drag/qS 
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lift coefficient, Lift/qS 


C L 

C m pitching-moment coefficient about c/2, Pitching moment/qSc 

Cp pressure coefficient 

D drag 

l airfoil maximum-thickness position, measured from leading edge 

L lift 

L/D lift-drag ratio, C L ^C D 

( L /D) m , v maximum lift-drag ratio 

M free-stream Mach number 

p free-stream static pressure 

q free-stream dynamic pressure 

S wing area 

t airfoil thickness 

t airfoil maximum thickness 

max 

t maximum thickness of portion of airfoil above chord line 

x projected length of airfoil surface on horizontal plane (parallel to free-stream 

flow) 

a angle of attack 

Subscripts: 

1 airfoil forward lower surface 

2 airfoil rearward lower surface 
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airfoil forward upper surface 


4 airfoil rearward upper surface 

MODELS AND APPARATUS 

Details of the four semispan wings are shown in figure 1. All wings were of delta 
planform with 65° leading-edge sweep and incorporated modified diamond airfoil sections 
having a thickness -chord ratio of 0.06. Airfoil geometry variables consisted of camber 
and longitudinal location of maximum thickness. Wing 1 employed the optimum modified 
diamond airfoil, based on two-dimensional shock-expansion theory, for a free-stream 
Mach number of 4. 63. The optimization procedure is discussed in the section "Theoreti- 
cal Considerations. " The wings were attached to a four -component strain-gage balance 
housed within a sting-mounted splitter plate. (See fig. 2. ) A gap of approximately 
0. 07 cm was maintained between the wing and plate to prevent fouling. 

Boundary -layer transition strips were placed 1 cm rearward (streamwise) of the 
wing leading edges. The strips were composed of No. 40 carborundum grains (average 
diameter approximately 0.046 cm) embedded in a plastic adhesive. The average diameter 
of the grains was approximately equal to the local boundary -layer thickness at M = 4. 63 
and exceeded the thickness by about 15 percent at M = 3. 95. Estimates based on the 
data of reference 8 indicate that the total particle -drag coefficients are within the accu- 
racy of the data. 

The investigation was conducted in the high Mach number test section of the Langley 
Unitary Plan wind tunnel, which is a variable-pressure, continuous -flow facility. The test 
section is approximately 1. 22 meters square by 2. 13 meters long. The nozzle leading to 
the test section is of the asymmetric sliding-block type, which permits a continuous varia- 
tion in Mach number from about 2. 3 to 4. 7. 

MEASUREMENTS, TEST CONDITIONS, AND ACCURACIES 

Aerodynamic forces and moments were measured at Mach numbers of 3.95 and 4. 63 
and at a Reynolds number of 9. 84 x 10 6 per meter. The angle-of-attack range was 
approximately -2° to 7°. The dewpoint was maintained below 239 K to minimize moisture 
condensation effects. 

The estimated accuracies, based on instrument calibration and data repeatability, 
are as follows: 
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± 0.0001 

± 0.002 

±0.0004 


THEORETICAL CONSIDERATIONS 


In the theoretical procedure employed herein, which utilizes the shock-expansion 

method, the flow is considered to be two-dimensional and hence affected only by airfoil 

shape. Figure 3(a) shows the effects of the longitudinal location of maximum thickness •£ 

t c 

on the maximum lift-drag ratio for various degrees of camber . u as predicted by the 

max 

two-dimensional shock-expansion method for a series of modified diamond airfoils having 
a thickness -chord ratio of 0. 06. The data are computed for inviscid flow and a free- 
stream Mach number of 4. 63. It should be noted that the peak values of maximum lift- 

drag ratio occur at — ~ 0. 60. A cross plot of these data, presented in figure 3(b), exhibits 
t c 

a peak value at t — — ~ 0. 60. Thus, the optimum modified diamond airfoil, in terms of 
t max 

maximum lift -drag ratio, has the maximum thickness located at approximately the 
60-percent-chord station and is cambered so that approximately 60 percent of the thick- 
ness is above the chord line. 


It was noticed during these calculations that for the more efficient airfoils, the pres- 
sure coefficients on the forward upper surface g and the rearward lower sur- 
face C 0 were very small at the inviscid (— \ ; for the optimum airfoil = 0. 60, 

. \ D/max \ c 

hi 


max 

zero. 


= 0. 60 , for all practical purposes, these pressure coefficients could be considered 


This trend is demonstrated in figure 4, which shows the corresponding variation of 
pressure coefficient and lift-drag ratio with angle of attack for the optimum profile (design 

M of 4.63, - ffy? = O.Oe). It is seen that at the maximum lift-drag ratio C 


Cp g are essentially zero. 


'P,3 


\ 

and 


This relationship indicates that at the optimum angle of attack, 
surface 3 is approximately parallel to the free-stream flow and, for the relatively thin air- 
foils considered, that surface 2 is also approximately parallel to the free stream. This 
relationship allows the use of a relatively simple optimization procedure instead of the 
laborious method employed to obtain the data of figure 3. In figure 5, surfaces 2 and 3 
are placed parallel to the free-stream flow. By assuming a wing- section width of unity, 
the following equation can be obtained: 
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c = ±- 

'"'L qS 


' V [( P l ' P ) X 1 + ( P 2 ' p ) x 2 - (»S - p )*3 - ( p 4 ’ P)^] 

_ Pi - P Xj f P2 “ P x 2 P3 “ P x 3 P 4 - p x 4 
qc q c q c q c 


*T *9 Xo x 4 

= C 1 — + C o — -C „ — - C 

p,l c p,2 c p,3 c p,4 c 


Since C p 3 = 0 and, for the relatively thin airfoils considered, C p 2 * 0 

X. X, 

C T = <T , -4 - C 


Thus 


p,l c ~ p,4 c 




cos a - x^ cos a ~ x z) 


- cos a — P — M — £i^(c cos a — 'j 

c cos a) c \ cos a) 


C L = C pjl ^cos « - L^j - C pj4 (cos « - 


(1) 


Similarly, 


C D = — 
D qS 


= ^ [( P l ' P ) ‘max cos “ - (»4 ' P) ‘max cos “] 

= W Msa .P4_2 W 

q C q c 


and 


'D “ ( ( 


Ct^ = c„ , - c 


t 


max 


P.l ~ ^4) c 


cos a 


( 2 ) 
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Since cos a is very nearly unity, equations (1) and (2), respectively, reduce to 


C 


L 



C 


D 




t 


max 


Hence 



(3) 

(4) 


(5) 


It should be pointed out that since surfaces 2 and 3 are assumed to be parallel to the free- 

stream flow, — and u are always equal. A plot of this relationship is shown in 
t c Wax 

figure 6 for m ^~ x = 0. 06 and free-stream Mach numbers of 4. 63, 6, 00, and 8. 00. It 

should be noted that for M = 4. 63, the peak value of lift-drag ratio occurs at ^ * 0. 60, 
which defines an airfoil of geometry essentially the same as the optimum shape obtained 
by the iterative procedure initially employed. For the higher Mach numbers the data of 
figure 6 indicate that the airfoil shapes are optimized at approximately equal to 0. 63 
and 0. 69 for Mach numbers of 6. 00 and 8. 00, respectively. Also of interest is the appre- 
ciable increase in the optimum lift-drag ratio with increasing Mach number. 


Figure 7 shows the variation of ( — | with 4 for Mach numbers of 6. 00 and 

\D/max c t 

8. 00, based on the shock-expansion method, for the optimum t values predicted by 

/ t i r max 

equation (5) is equal to - if surfaces 2 and 3 are assumed to be parallel (see 

\ Vmax c 

fig. 5)1. Calculations for M = 4.63 are included for comparison. The symbols repre- 
sent the optimum configurations according to equation (5). Similar data are presented 

in figure 8 for the variation of (^] with , tu , based on the shock-expansion 

w/max tmax 

method, for the optimum ^ values as predicted by equation (5). Although analytical 
data were calculated for limited ranges of maximum-thickness location and camber 
at Mach numbers of 6. 00 and 8.00, the data of figures 7 and 8 indicate that equation (5) 
also provides close agreement with the shock-expansion method in predicting the opti- 
mum modified diamond airfoil at these higher Mach numbers. 
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COMPARISON OF THEORY AND EXPERIMENT 

Comparisons of the two-dimensional shock-expansion predictions with experimen- 
tal data at Mach numbers of 3.95 and 4. 63 are presented in figures 9 to 12 for wings 1 
to 4, respectively. It should be emphasized that the theoretical method considers the flow 
to be two dimensional and hence affected only by airfoil shape. In all cases the predicted 
pitching- moment coefficients are in excellent agreement with the experimental values. 
Theory underestimates lift coefficient, with the difference increasing with angle of attack. 
At angles near cruise attitude (3° « a « 5°) the decrement ranges from 3 to 7 percent. 

The predicted drag coefficients are considerably less than the experimental values, 
with the decrement in minimum drag coefficient ranging from about 12 to 17 percent at 
M = 4.63 and from about 15 to 19 percent at M = 3.95. The theoretical values include 
viscous effects based on the T’ method of reference 9 and the assumption that fully devel- 
oped turbulent flow exists over that portion of the wing rearward of the boundary -layer 
transition strips. The results of the evaluations reported in references 10 and 11 indi- 
cate that the theoretical skin-friction coefficients predicted by the T’ method of refer- 
ence 9 agree well with experimental data in the Mach number range considered herein. 

It is believed that the underestimation of drag coefficient is due primarily to the 
unaccounted-for contribution of the three-dimensional effects of wave drag. 

A summary of the theoretical and experimental maximum lift-drag ratio is presented 
in figure 13. Theory provides an accurate estimate of the relative effects of maximum- 
thickness position and camber on maximum lift-drag ratio; however, maximum lift-drag 
ratio is overestimated by approximately 6 to 8 percent. Although the amount of experi- 
mental data is limited, the results of figure 13 suggest that for the Mach numbers of this 

investigation and - — — = 0.60, the optimum location of maximum thickness is in the 60- 
t max 

to 65 -per cent- chord range and agrees closely with the theoretical predictions. Insufficient 
test data are available to indicate the optimum camber, but the values obtained exhibit 
trends paralleling those of theory. Hence, although additional data might be desired to 
define the optimum airfoil geometry more clearly, the data obtained show that the two- 
dimensional shock- expans ion method is very useful in estimating the optimum airfoil shape 
for a finite wing at high supersonic and hypersonic speeds. 

CONCLUSIONS 

An investigation has been conducted to explore the potential for optimizing airfoil 
shape at high supersonic speeds by utilizing the two-dimensional shock-expansion method. 
Theoretical data are compared with the experimental force and moment coefficients 
obtained on four delta-planform semispan wings having a leading-edge sweep angle of 65° 
and incorporating modified diamond airfoils with a thickness -chord ratio of 0.06. The 
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wings differ only in airfoil maximum-thickness position and camber. The experimental 
data are obtained at Mach numbers of 3. 95 and 4. 63 and at a Reynolds number of 
9.84 x 10® per meter. The conclusions are as follows: 

1. A relatively simple method is developed for predicting, in terms of lift-drag 
ratio, the optimum modified diamond airfoil at high supersonic and hypersonic speeds. 

2. Although maximum lift-drag ratio is overestimated by approximately 6 to 8 per- 
cent, theory provides an accurate estimate of the relative effects of airfoil maximum- 
thickness position and camber. 

3. Theory underestimates lift coefficient, with the difference increasing with angle 
of attack. At angles of attack near cruise attitude (3° * a ~ 5°) the decrement ranges 
from 3 to 7 percent. 

4. The predicted drag coefficients are considerably less than the experimental 
values, with the decrement in minimum drag coefficient ranging from about 12 to 17 per- 
cent at a Mach number of 4. 63 and from about 15 to 19 percent at a Mach number of 3. 95. 

5. The predicted pitching-moment coefficients are in excellent agreement with the 
experimental values. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va. , April 3, 1974. 
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l/c=0. 50 l/c=0. 60 l/c=0. 70 


Figure 1.- Model details. Linear dimensions are in centimeters. 




Figure 2.- Tunnel installation 
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(a) Effects of If c at constant tjt max . 

Figure 3.- Maximum lift-drag ratio as predicted by two-dimensional shock -expans ion method. Inviscid flow; 

M = 4. 63; t/c = 0, 06. 
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(b) Effects of t u /t: max at constant Z/c. 
Figure 3.- Concluded. 
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Figure 4.- Variation of pressure coefficient and lift-drag ratio with angle of attack for 


optimum profile. Design Mach number of 4. 63; inviscid flow; = 0. 06 
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Figure 6.- Variation of L/D with Z/c, based on equation (5). 
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Inviscid flow; 
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= 0.06. 
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Figure 8.- Variation of (L/D) max with t u /t max , based on shock-expansion method, 
for optimum l/c values predicted by equation (5). Symbols represent optimum 
airfoils according to equation (5). 
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(b) M = 4.63. Concluded. 
Figure 11. - Concluded. 
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Figure 13,- Summary of experimental data and two-dimensional shock- expans ion predictions. 







